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abstract 

We study a generality of an inflationary scenario by integrat- 
ing the Einstein equations numerically in a plane-symmetric 
spacetime. We consider the inhomogeneous spacetimes due to 

(i) localized gravitational waves with a positive cosmological 
constant A, and (ii) an inhomogeneous inflaton field $ with 
a potential ^m 2 ^ 2 . For the case (i), we find that any ini- 
tial inhomogeneities are smoothed out even if waves collide, 
so that we conclude that inhomogeneity due to gravitational 
waves do not prevent the onset of inflation. As for the case 

(ii) , if the mean value of the inflaton field is initially as large 
as the condition in an isotropic and homogeneous inflationary 
model (i.e., the mean value is larger than several times Planck 
mass), the field is soon homogenized and the universe always 
evolves into de Sitter spacetime. These support the cosmic no 
hair conjecture in a planar universe. We also discuss the effects 
of an additional masslcss scalar field, which is introduced to 
set initial data in usual analysis. 
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field under the existence of initial thermal fluctuations using new inflation- 
ary model. They found that inflation occurs only when the potential is flat 
enough. Since their motivation was different from studying the cosmic no 
hair conjecture, they did not consider the case which include large inho- 
mogeneities of spacetime. Although they found de Sitter space is a final 
state, their initial inhomogeneity is not enough large to study the cosmic 
no hair conjecture. 

Assuming a spherically symmetric spacetime, Goldwirth and Piran (GP) 
studied the behavior of inhomogeneous distributions of scalar field. From 
an estimation of an energy balance between potential energy and scalar 
field gradient term, they proposed a criterion for the onset of inflation 
such that 

a scale length of inhomogeneity 5<fi 

- — : > a few times (1) 

horizon scale m p i 

where 5<fi is a spatial deviation of the scalar field from the mean value 
and rripi is Planck mass. They found that a short scale inhomogeneity 
compared to the horizon scale leads spacetime to collapse in closed uni- 
verse, and concluded that is sufficient condition for chaotic inflationary 
models. GP's simulations, however, are done under assumptions on the 
isocurvature initial data with additional massless scalar field, which may 
cause homogeneous expansion of the spacetime in the very initial stage. 
This additional scalar field may change a homogenization process of space- 
time as we will see later. 

In this paper, we study the validity of the cosmic no hair conjecture in a 
plane symmetric spacetime. We have mainly two reasons why we assume a 
planar universe. First, we investigate the behaviors of the inhomogencitics 
due to the gravitational waves, which cannot exist in a spherically symmet- 
ric spacetime. It is well studied that any linear gravitational waves decay 
and disappear immediately in de Sitter spacetime so our main inter- 
ests are the spacetimes with strong gravitational waves. As some numerical 
study in flat spacetimes shows, if the strong gravitational waves are local- 
ized, those come to have apparent horizons and collapse into black holes 



We show the effects of inhomogeneities due to gravitational waves in Sec- 
tion 3, and those due to inflaton field in Section 4. Conclusions and remarks 
are in Section 5. As units, wc use 8nG = 1 and c = 1. 

2 Equations and Numerical Methods 
2.1 The models 

We consider following two effects on the evolution of the spacetime; [I] 
inhomogeneities due to gravitational waves and [II] inhomogeneities due to 
an inflaton field. To make the problem simple, we set models as follows: 

[I] Initial large inhomogeneities due to gravitational waves exist in the 
spacetime with positive cosmological constant A (no inflaton field). 
We are interested in the competitive processes between the expansion 
of the background spacetime and the attractive forces of the gravity. 

[II] Introduce an inhomogeneous inflaton field <& in the context of the 
chaotic inflationary model. In a realistic inflationary model, it is 
natural to assume such an inhomogeneous inflaton field. Our interest 
is the effects of the spatial gradient term of the field on the evolution 
of the spacetime, which may prevent the onset of the inflation. 

For the second model, we make initial data in two ways. One is the 
similar setting of initial data as those of GP's, i.e., the spacetime with an- 
other massless scalar field to obtain an isocurvature initially. The other 
is without such an additional scalar field. We examine the contributions 
of \I/ to evolutions of the spacetime. 

In the chaotic inflationary scenario, we usually assume that the average 
of the inflaton, $, is larger than a few m p i in order for the Universe to 
expand sufficiently. Even under an existence of large inflaton field as an 
average ($ > a few m p i), it is not trivial how the large inhomogencity 
of <& contributes to the evolution of the Universe. We analyzed whether 
or not such an gradient term in the energy p$/ prevents from the onset of 
inflation. 



where h v — 6" + n u n^ is the projection operator onto E. 

In section 3, we will treat vacuum spacetime, i.e. T^ v — 0, under the 
assumption of the existence of a positive cosmological constant A, while in 
section 4, we will treat a scalar field <3? with potential V($) = ^m 2 $ 2 as 
a source of the inflation rather than A. The energy momentum tensor of 
the scalar field is given by 

Tpu = ft,*ft* - 5 ^(V<f) 2 + V{$)]. (10) 
The scalar field evolution, then, obeys the Klein-Gordon equation 

To integrate (|ll]), we introduce the conjugate momentum 

n = ^(-a t $ + -?-d x $), (12) 

a 7n 

where 7 = det 7^ , and write down eq. ( |TT| ) into two first-order partial dif- 
ferential equations: 

ft* = — ft*--^n, (13) 

7ii \/7 

9 t n = a^7^ + ft — [^n-aV7^*]. (14) 
a<P 711 

The dynamical variables are 7^ and Ify (and $ and n, when a scalar field 
exists). 

We impose periodic boundary condition, i.e., f{x\,t) = /(xn+i,^), 
where N is the grid number in x-direction. 

Our simulation procedure follows Nakamura, Maeda, Miyama and Sasaki 



(i) Determine initial values by solving two constraint equations. 

(ii) Evolve time slices by using the dynamical equations. 

(iii) Check the results of (i) and (ii) using the two constraint equations on 
every time slice. 



In order to get the periodic solution of (|l8|), the value of trK which 
we give on the initial hypersurface should be chosen so that such a so- 
lution exists. The condition for trK depends on the models. Details are 
shown in the following sections. We use the 5th order Runge-Kutta method 
(Fehlberg method) to solve (j) in (|l8|). 

2.4 Time Developing Scheme 

For coordinate conditions, we impose a geodesic slicing condition such as 
a = 1 and /3 = 0. This slicing condition entails the risk that our numerical 
hypersurface may hit a singularity and stop there. If no singularity appears, 
however, then it may be the best coordinate condition for revealing whether 
de Sitter space emerges as a result of time evolution. 

We use finite differential scheme to integrate the Einstein equations (0) 
and (||) with 400 grids in the spatial direction x. We use a simple central 
finite difference with second order scheme to determine derivatives. 

The time step At on each slice is determined, as to keep the accuracy of 
the computations. We set At by imposing that a volume element 7 and a 
expansion rate trK do not change too much (±0.1% change) in each step. 

The accuracy is checked by comparing both sides of constraint equations. 
For Hamiltonian constraint equation in the case of vacuum, for example, 
we calculate 

V ' ' \R^)\+K 2 + K l] K*J + 2A+|2p H | V ; 

at the point x on £(t) and regard its maximum value Err or (t) — 
iriax{Err(t, x)\ x € £(£)} on each hypersurface as the error on that slice. 
A similar definition is performed for momentum constraint equation. An 
example of Error(t) is shown in Fig. 2.1. In our calculation, the maximum 
error in the Hamiltonian constraint equation is a few times of O(10~ 4 ) on 
the initial hypersurface, and this accuracy is maintained even after time 
evolutions. 

Our overall procedures to evolve the system are as follows. At first stage, 



As we already mentioned in the introduction, linear gravitational waves 
are always decay in the de Sitter background spacetime. This result is 
also confirmed in our numerical calculations as one of our code tests. We, 
therefore, concentrate our attentions to the large inhomogeneities due to 
gravitational waves. The summary has been reported in p3| . 

Since we have A in this system, we adopt the Hubble expansion time 
th = (A/3) -1 / 2 as our time unit, which is a characteristic expansion time 
of the universe. Our unit of length is also normalized to the horizon length 
of the de Sitter spacetime Ih = (A/3) -1 / 2 . 

To evaluate the inhomogeneities on each hypersurfaces, we use three 
different curvature invariants. First one is the 3-dimensional Riemann 
invariant scalar ( 3 'Rijkl ^ 3 'R l ^ kl , where ( 3 'Rijkl is the Riemann tensor 
of the 3-metric on S. We use its dimensionless value normalized by the 
cosmological constant, 



C(t 7 x) — — ±— on L(i), (21) 

and call it the "curvature" hereafter. We estimate the magnitude of the 
inhomogeneities in the 3-space £ by the maximum value of this "curvature" 
on each slice, i.e. C max (t) = max{C(i, x) \ x e £(£)}. 

Second one is an invariant scalar induced from the Weyl tensor C^ vp a- 
For the vacuum spacetime, if the Weyl tensor vanishes and no singularity 
appears, the spacetime is homogeneous and isotropic. Hence this can be 
used to check whether de Sitter universe is recovered or not. We use the 
decomposition of the Weyl tensor: 

Ep a C p pal ,Tl^Tl , B pa — C ppavTl^Ti , (^) 

where *C ilvpa = jC'^^Sap^ is the dual of the Weyl tensor and is 
a timelike vector orthogonal to the hypersurface S. In analogy to the 
electromagnetism, the 3-dimensional variables E pa and B pa are called an 
electric and a magnetic parts of the Weyl tensor. We can reconstruct the 
Weyl tensor completely from this pair of tensors. Moreover, those two 



on the initial hypersurface T,(t — 0). Here, 8k is a positive constant, which 
should be introduced because the periodic pulse increases expansion rate 
of the universe besides A. By integrating ([l8]), we get conformal factor cf>. 
There is one free parameter A, which fixes the scale. If <fi is multiplied by A, 
we find another solution by replacing 8 k with 8k /A 2 , which is physically 
the same as the previous solution. Hence, if we impose — > 1 at the 
numerical boundary, then the value of 8k is fixed; it is less than O(10 -2 ). 

A pulse-like wave has two characteristic physical dimensions, a width 
and an amplitude. We use C max (i) as an amplitude measure, and we define 
the width l(t) by the proper distance between two points where 7 (the 
square of the 3- volume) decreases by half from its maximum value ^max- 
An example of the initial data is shown in Fig. 3.1. This is the case of 
l(t = 0) = 0.086Z H ,C max (f = 0) = 9.67. 

The behaviors of the time evolutions are following. Initial distortion 
begins to propagate both in the ±z directions. We see the initial curvature 
C(t = 0,x) localized at the center with three peaks [Fig. 3.1(b)] separates 
just after t = 0, and propagate at the light speed. In Fig. 3.2, we show 
both "curvature" C(t,x) and "super-energy" H(t,x). Two waves moves to 
the numerical boundaries and collides each other on the boundary, since we 
assume periodic boundary condition. After the collision, two waves proceed 
again in their original directions. This behavior is reminiscent of solitonic 
waves. We also show, in Fig. 3.3, C max (i) and the invariant X(t, x) = 
V 1 RiivpoR^ 90 ' I 'A. We can see clearly that the collision of the waves occurs 
around t = 0.2th, when the "curvature" is superposed. At the final stage, 
the spacetime is homogenized by the expansion of the universe. We can 
see T(t, x) — > y/8/3 (homogeneous de Sitter spacetime) within one Hubble 
expansion time th as expected [see eq.(p6|)]. We see the spacetime finally 
succumbs to the overall expansion driven by the cosmological constant, 
and becomes indistinguishable from de Sitter spacetime. 

In our simulation, I and C max (0) range for [case 1] between 0.080/if < 
I < 2.51h and 0.020 < C max (0) < 125.0. We also computed the evolution of 
[case 2] , and found that all initial inhomogcncitics decay and disappeared 



lutions for the case of colliding waves and for the case in a contracting 
background are still under studying. The details will be reported in the 
elsewhere. 

4 Inhomogeneities due to scalar field 

In this section, we introduce a scalar field &(x,t) as a source of inflation 
instead of cosmological constant, and show the results of the behavior 
of its inhomogeneities to the evolution of the spacetime. In the usual 
analytic approach of inflationary scenario, the infiaton field is assumed 
homogeneous. By introducing inhomogeneities, the energy density of the 
scalar field /?$ is written as 

p* = ^n 2 + -^(^$) 2 + V ($) 
27 2711 

= Pi + pv+V^), (30) 

where V($>) is the potential term of the scalar field. Our interests is how 
the gradient term p$i affects on the expansion of the universe. The initial 
kinetic term p^ gives less effects on the inflationary scenario. 

The same-aimed simulations have been already done by GP using a 
spherically symmetric code 0, and one of their results is that the space- 
time with large p& do not evolve into the inflationary era. In order to 
compare the results between ours and theirs, i.e., the difference of the 
symmetry of the spacetime, we first set similar initial situations and pa- 
rameters with GP, that is with 'isocurvature initial data' (section 4.1). 
This condition makes us easy to determine initial data with constant mean 
curvature slicing, but also makes us unclear the effects of infiaton <f> itself, 
because the spacetime may expand homogeneously in the very initial stage 
by introducing another scalar field and such an expansion might reduce 
an initial large inhomogeneity of infiaton field. So that we also prepare the 
initial data without and computated their time evolutions (section 4.2). 



We choose $ several times of m p \ in order to get 'suitable' expansion of 
the Universe when = 0. This is from a condition for inflation in isotropic 
and homogeneous universe. We are interested in effects of gradient term 
of the scalar field p$' to the time evolution of spacetime, especially large 
compare to ptotai- We get large p$<, if we choose small A and large SQ. 
However, the total energy density ptotai = P$' + V{<f>) at the inflationary 
era is order of the Planck scale m^, so that if we introduce 'suitable' <E>o, 
then acceptable values of A and 5$ are limited. But we can get large 
enough p$/ compare to ptotai as will be shown in Table. 

First, we prepare initial data of the scale A ~ Ir and the amplitude up 
to <5<I> ~ 0.1$o (model [Il-la ~ lc] in Table 4.1). Corresponding ratio of 
Pv = +P*' to ptotai ar e also shown in Table 4.1. Those model have large 
contributions of pv, we may regard that gradient term is locally dominated. 
As we can see from Fig. 4.1, in which the ways of evolutions (model [II- lb]) 
is shown, all the inhomogeneous scalar field are first homogenized and the 
field go to the slow-rolling phase along to the ordinary chaotic inflationary 
scenario. We find the spacetime expands sufficiently in the region where 
the inflaton has large enough <!>o. This confirms the criterion ([!]) proposed 
in spherically symmetric case by GP. 

Next, we prepare initial data of the small A. In these cases, we may 
expect a formation of the collapse of the spacetime, as an example shown 
by GP. Our case has planar symmetry, so if the collapse occurs, then the 
naked singularity appears. Hence from this calculation we can also discuss 
the so-called cosmic censorship in the spacetime with inflaton field. We 
prepare data with the same 8& and different A's such as A = Ih, ^Ih, \Ih 
and gZff, when $o = 8.0m p i and 5$ — 0.06m p i (model [Il-ld ~ lg] in 
Table 4.1). In the last case (A = gijj), Pv is locally 86% of the p to tab we 
can regard gradient term is locally dominated again. 

The typical time evolution of $ is shown in Fig. 4.2 (model [Il-lg] ) . The 
evolution is quite similar to Fig. 4.1. That is, in all the case, the field $ is 
homogenized and the spacetime expands sufficiently. 

For various scale of the deviation A (model [Il-ld ~ lg], we show the 



4.2 Evolutions of the inhomogeneous scalar field 

Next, we show the time evolution of the initial data without introducing a 
massless scalar held. This setting makes clear the effects of inhomogeneous 
scalar field, since initial expansion rate are not uniform over the spacetime. 

To show the difference with the previous model, we set the same chaotic 
potential with m 2 = 0.01, and the same initial scalar distribution form 
^initial- To define a unit, we set the effective cosmological constant as 
A ff = im 2 $oi where &q is the mean value of the initial scalar field [see eq. 
(p2|)], and fix unit Ih and tu using A e g. We set the background spacetime 
is flat on the initial hypersurface, and choose tvK as 



Since our present model is not isocurvature, 8k must be introduced in 
order to take into account the inhomogeneity effect, and is determined by 
the same condition as in Sec. 3. 

As like the cases in the previous section, we first prepare initial data for 
large amplitude <5<£> and the scale A to be larger than a horizon scale (model 
[II-2a ~ 2c] in Table 4.2). Initial rate of p$/ to ptotai and each 8k are shown 
in Table 4.2. The gradient energy rate pi^/ptotai are large enough to treat 
our problem. The inflaton in such a case behaves like a homogeneous field 
in one horizon region, and induce a sufficient expansion if the average value 
of $ is 'suitably' large for the inflationary model. These results are again 
agree with GP's results as we expected. 

We next prepare initial data for small A compare to the horizon scale. 
With the same reason in the previous section, if we concentrate our atten- 
tions for the small scale deviations A compare to the horizon scale, then 
only small amplitude of 8$ is acceptable. We set A = 0.2i# and change 
8<S> (model [H-2d - 2h] in Table 4.2). 

A typical example of the time evolution of scalar field is shown in Fig. 
4.7. We see the field $ is rolling down along the chaotic inflationary po- 
tential. In the figure, only the configuration until 12t# is drawn, but we 
can see differences with that of the isocurvature model (Fig. 4.2). Since 




(34) 



time; strong inhomogeneity [I] due to localized gravitational waves in the 
spacetime with positive cosmological constant and [II] due to inhomoge- 
neous inflaton field. 

For the model [I], we see all initial inhomogeneities are smoothed out 
within one Hubble expansion time tjj including their collisions' effect. So 
that we conclude that the nonlinearity of the gravity do not cause to pre- 
vent the onset of the inflation under the existence of positive cosmological 
constant A. 

For the model [II] , we prepare two initial situations; isocurvature initial 
data and fully inhomogeneous initial data. In both cases, we find that 
initial gradient term of the energy density of the scalar field may not lead to 
collapse of the spacetime. If the initial mean value of the inflaton field, $o, 
is large enough to derive inflation in homogeneous and isotropic spacetime, 
then we can conclude such a universe will always be homogenized and 
expand sufficiently. 

For the case of the inhomogeneous scale A is larger than the horizon 
scale, our calculations confirm the criterion ([!]) by Goldwirth and Piran, 
and also their calculations for spherical spacetime. 

Even for the case of small A, we could not find the collapsing universe. 
We did not face to a formation of a naked singularity at least in a planar 
universe. This may support a cosmic censorship conjecture. Therefore, 
we conclude that inflation is generic, and the usual analysis with isotropy 
and homogeneity may be justified in a plane-symmetric spacetime. We 
also partially confirm our inflationary scenario of inhomogeneous universe, 
which we mentioned in the introduction. That is, a large inhomogene- 
ity does not necessarily prevent cosmological constant-driven expansion, 
instead the spacetime evolve into a de Sitter spacetime with many small 
black holes induced by gravitational collapses. 

As we denoted in section 4.2, the condition of isocurvature initial data 
by introducing an additional massless field may lead to the different con- 
figuration on the evolution of the spacetime. The critical difference did 
not appear in the plane symmetric spacetime, but it is necessary to treat 
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Table Captions 



Table 3.1: 

Examples of initial parameters for model [I] [case 1]. a and xq are param- 
eters to define jij in eq.(p7f), L is a periodic boundary length, and 8 k is 
a constant in eq. (p9|) . l(t — 0) and C max (t = 0) are width and amplitude, 
respectively, on the initial hypersurface. 

Table 3.2: 

Examples of initial parameters for model [I] [case 2]. a and xq are pa- 
rameters to define jij in cq. (p8|) , and the rest are the same with Table 
3.1. 

Table 4.1: 

Examples of initial parameters for the isocurvature model. $0, S^>, A means 
average, amplitude of deviations and scale of the deviation, respectively [in 
cq.(|2^)]. Initial ratio of gradient term pv to the ptotai are also shown. The 
mass term of the potential, m, is set as m 2 — 0.01. 

Table 4.2: 

Initial parameters for the fully inhomogeneous model. 5k is a constant in 
eq.(p4). The mass term of the potential is set as m 2 — 0.01. 



Model 


$ 


(5$ 


A 


(max) Pv/Aotai 


Pv/ Ptotal 


Model 


$ 


6$ 


Il-la 


5.00 


0.75 


1.00 


92.5% 


46.1% 


II-2a 


8.00 


1.50 


II- lb 


8.00 


0.72 


1.00 


85.3% 


42.5% 


11-26 


8.00 


2.40 


II-lc 


8.00 


0.36 


0.50 


85.3% 


42.5% 


II-2c 


8.00 


3.60 


Il-ld 


8.00 


0.06 


1.00 


2.37% 


1.18% 


II-2d 


8.00 


0.03 


Il-le 


8.00 


0.06 


0.50 


9.47% 


4.73% 


II-2e 


8.00 


0.06 


Il-lf 


8.00 


0.06 


0.250 


37.9% 


18.9% 


11-2/ 


8.00 


0.18 


Il-lg 


8.00 


0.06 


0.167 


85.2% 


42.5% 


II-2 3 


8.00 


0.30 


II-2/i 


8.00 


0.60 



Tabl. 

Table 4.1 
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Figure Captions 



Figure 2.1: 

The accuracy of the Hamiltonian constraint equation. The solid line and 
the dotted line show the maximum error Error{t) and the averaged error, 
respectively, on the hypersurface at every time step for the calculation of 
the inhomogeneous scalar field (the case of model [II-2f] in Section 4.2). 

Figure 3.1: 

An example of the initial configuration for the case of a pulse-like distortion 
([case 1]). This is the case of model [I-la] in Table 3.1. In (a), the conformal 
factor (f> in equation (|l5|) is shown. In (b), the "curvature" C(t,x) in 

Figure 3.2: 

Time evolutions of propagating plane waves, (a) "curvature" C(t,x) 
[cq. (pi])] and (b) "super-energy" Ti(t,x) [eq.(p5|)] are shown for the ini- 
tial data shown in Fig. 3.1. 

Figure 3.3: 

The maximal value of T(t, x) = R^ vpa R^ v P' J / K (solid line), and C(t, x) = 
y/ ^Rijki ( 3 )ii* jfei /A (dotted line) on each hypersurface £(i) are shown for 
the same data with Fig. 3.1. We find 2(t,x) — > a/8/3 within one Hubble 
expansion time (homogeneous de Sitter spacetime) as expected. 

Figure 3.4: 

The time evolution of the "curvature" C(t, x) resulting from waves that are 
located closely. Two waves are the same form, / = 0.10Z# and C max (£ = 
0) = 51.0, and the periodic distance is 0.30Z#. We see them collide and 
in the collision region, the "curvature" seems to be superposed, but finally 
the spacetime is homogenized by the expansion of the universe. 



Figure 4.1: 



The proportions of the each terms in energy density to ptotai I /°$ (bold-solid 
line), (bold-dotted line), p^, (thin-solid line), p^i (thin-dotted line), and 
y($) (three-dot-dashed line), (a) shows until 1.0th, while (b) shows until 
5.0th- 

Figure 4.6: 

The evolutions of the maximum Weyl invariant ^C^paC^P 17 / A e g on the 
hypersurface. 

Figure 4.7: 

A typical example of the time evolution of initially inhomogeneous scalar 
field. The model [II-2f] in Table 4.2 is shown. We see the field <E> is rolling 
down along the chaotic inflationary potential V^($) = ^m 2 ^ 2 with m 2 = 
0.01. Compare to the Fig. 4.2, we notice the scalar field's configuration 
drags its initial form even after the time evolution. 

Figure 4.8: 

Trajectories in $ — II phase space for various kinetic terms' initial data for 
the fully inhomogeneous model (at x — 0, where the maximum $ exists 
initially). Only the beginning part is drawn. The detail parameters are 
shown in Table 4.2. The bold line shows the trajectory in the case of 
homogeneous scalar field ($ = 8.0m p ;). The dotted line, dot-dashed line, 
dashed line and solid line are those of Model [II-2e, 2f, 2g, 2h] respectively. 
All of them will coincide to the homogeneous' line soon after slowing down 
the potential. 
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